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Part  I  -  An  Approximate  Treatment  (S.  V.  Yadavalli) 


SUMMARY 


V 


Assuming  a  decelerated  charged  particle  beam  penetrating  a  uniform 
plasma,  a  relation  for  the  growth  rate  of  electrostatic  instabilities  is 
obtained.  It  is  shown  that  the  growth  rate  in  a  decelerated  beam  case 
is  always  higher  than  in  a  uniform  one.  An  experimental  method  of  de¬ 
termining  the  ^electron^ plasma  frequency  of  the  plasma  and  the  deceler- 
/■« I4*-  \lj.h  2 

ating  parameters^  and  (J  is  also  suggested. 
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I. 


Introduction 


In  this  report  we  treat  the  problem  of  the  interaction  of  a  charged 
particle  beam  with  a  plasma,  where  the  plasma  was  created  originally  by 
the  passage  of  the  charged  particle  beam  through  a  neutral  gas.  This 
problem  is  of  importance  in  connection  with  the  experimental  work  now  in 
progress  at  the  University  of  Illinois.  Also,  mention  must  be  made  that 
some  experimental  work  along  similar  lines  (involving  the  stability  of 
charged  particle  beams)  is  being  carried  on  at  the  Research  Laboratory 
of  Electronics  at  M. I.T.  under  the  direction  of  L.  D.  Smullin. 

We  will  study  only  one- dimensional  situations  here  with  the  hope 
that  we  might  be  able  to  extend  the  treatment  later  to  include  the  case 
of  a  finite  beam  with  axial  deceleration. 

II.  Theoretical  "Model" 

Since  in  our  problem  the  charged  particle  beam  creates  its  own  plasma, 
one  expects  the  energy  of  the  beam  to  decrease  along,  say,  the  direction 
of  motion  of  the  beam.  Such  a  decrease  in  the  energy  of  the  beam  is  in¬ 
corporated  here  according  to  two  models,  by 
Case  I: 

assuming  that  the  charged  particle  beam  is  decelerated  exponentially, 
that  is,  in  the  non-relativistic  regime 

u(z)  ~  u  e  02  (1) 

o 

where  u  is  the  (entrance)  velocity  of  the  beam  at  z  =  0, 
o  ' 

a  is  some  positive  constant  determinable  by  experiment,  and 
z  is  the  pertinent  spatial  coordinate  along  which  deceleration 
is  experienced. 
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We  assume  further  that  the  steady  beam  current  is  invariant  with  z. 
Also,  the  plasma  thus  created  is  assumed  to  be  uniform,  and  the  plasma 
properties  are  independent  of  z. 


Case  II: 

assuming  that  the  charged  particle  beam  is  decelerated  according  to 


u(z)  ~  uq 


(i  *  ? 


(2) 


where  f3  >  0,  and  a  >  0. 

Note  that  according  to  the  model  of  Case  II,  the  beam  has  some 
finite  energy  at  z  =  oo,  because 


u(z) - »  uq(1  +  (3),  and 

lim  z  — »  0 

u(z) - )  uq 

lim  z  — >  » 


It  will  be  further  assumed  in  this  case,  just  as  in  Case  I,  that  the 
beam  current  is  invariant  with  z,  and  that  the  plasma  properties  are 
uniform  and  independent  of  z, 

III.  An  Approximate  Treatment 

In  this  part  of  the  report  we  will  consider  only  Case  II.  In  Part  II 
of  this  report,  H.  H.  C.  Chang  has  treated  both  the  cases  in  a  more 
rigorous  fashion. 

Let  us  now  proceed  with  the  discussion  of  Case  II.  We  assume  a  de¬ 
celerating  (non-relativistic)  electron  beam  penetrating  a  plasma.  We 
assume  here  that  the  charged  particle  beam  is  of  such  low  density  that 


All  transverse  variations  of  quantities  will  be  neglected  here. 
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it  does  not  affect  the  properties  of  the  "plasma"  medium  through  which  it 
traverses.  It  is  assumed  that  the  beam  and  the  plasma  are  "well  mixed." 
Let  the  velocity  of  the  beam  be  u;  then, 

u(z)  =  UQ  [^1  -t-  p  e  °®j  +  vb(z,  t)  (3) 

where  v(z,  t)  is  the  perturbation  in  velocity.  Assuming  only  electric 
fields,  we  may  write  the  force  equation  in  linearized  form  as 

3T  +  [uol 1  +  p  e_QZ)]  (  v  '  *b  +  v  ■  (uoP  e_Q2])  "  -’»*  (4) 

Assuming  that  all  perturbation  quantities  vary  as  eiwt  ikz 
Eq.  (4)  may  be  written  in  the  form 


iuivb  +  UQJl  +  p  e  02 1 1 -ikvb  -  uoC«p  e  02  j  =  -  tj[E(z)  +  E(z,  t)]  (5) 

where  E(z)  is  the  "equivalent"  electric  field  causing  the  deceleration, 
and  £(z,  t)  is  the  perturbation  in  the  electric  field. 

Then, 

iwvb  +  uq(  1  +  p  e  |-ikvb)  =  -  t]E  (5a) 

and  I  i  I  i 

u  1  +  3  e  02  -.u  op  e  ^  =  -  r|E(z)  (5b) 

o'  *  o  f 

Eq.  (5a)  is  the  one  of  interest  to  us,  which  may  also  be  written  as 


0  =  _ ini _ 

b  ju  -  kUol1  +  P  e_QZ)} 

A 

Using  the  equation  of  conservation  of  charge,  V  •  Jb 


(6) 

^b 

St  ’ 
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and  the  fact  that  the  unperturbed  beam  current  is  independent  of  z, 

★  l  —QZl  “  i 

which  means  that  p(z)  =  p  11  +  p  e  ,  we  obtain 


or 


_  "o  br  T  K  w  i 
b  -  kuQ  1  +  p  e  °*J  1 


A  i 

— QZ 

“  X 

-  iriE  k  p 

O  ' 

l  +  M 

<  u  -  ku 

i  o1 

,  „  -QZ ’ 

[72 

1  +  P  e 

11 

(7a) 


(7b) 


Eq.  (7b)  is  obtained  by  substituting  Eq.  (6)  in  Eq.  (7a). 

A 

dimensional  treatment  here  it  is  legitimate  to  write  E  = 
and  obtain 


k2  n^Pol1  +  p  e_QZ| 

[w  -  kuo[l  +  p  e-QZ)]2 


In  the  one- 

A 


** 

Sz 


(8) 


The  following  comments  are  in  order.  In  the  plasma  we  assume 
equal  numbers  of  ions  and  electrons,  and  as  the  ions  are  considerably 
heavier  than  the  electrons,  we  assume  that  only  the  electrons  in  the 
plasma  contribute  to  any  current  due  to  their  motion.  Or, 


(p)  A 
p  v 

°  P 


A 


(9) 


where  PQ^P^  is  the  unperturbed  charge  density  (of  electrons  alone)  in 
the  plasma.*** 

* 

p(z)  is  the  charge  density  of  the  beam. 

**  jr 

Cp  is  essentially  an  electrostatic  potential. 

*** 

This  charge,  however,  is  assumed  to  be  neutralized  in  the  plasma  by 
ions.  A. so,  Vp,  p^,  and  J  are  the  perturbed  velocity,  electron 
charge  density  and  current  in  the  plasma  respectively. 
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Employing  the  equation  of  continuity 


ikp 

o 


(P)  A 

P 


we  obtain 


From  the  force  equation,  we  have 


or, 


One  may  then  write  Poisson's  equation 


(10) 


(11) 


(12) 


^  <J)  =  +  ’  a,ld  °btaitl 


after  some  simple  algebra  that 


rip  /€ 

1  O  0 


11  +  g  e 


~az 


-l 


[a)  -  kujl  +  p  e_0!Z)] 
After  a  little  algebra,  we  find  that 


k  = 


03 

»  o 


% 


HP, 


(P) 


e  co 
o 


,3/2  ,,  2,2. 

U  t  (1-03  /03  ) 

O*  P 


(13) 


(14) 


where  £ 

= 

1,  „  -az\ 

[l  +  P  e  j 

2 

% 

= 

_,1po/eo 

2 

0) 

p 

= 

-Tjp  ^pVe 

1  o  o 

and 

CDp  are  the  "electron"  plasma  frequencies  of  the  beam  and 

the  plasma,  respectively. 
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Eq.  (14)  is  the  desired  relation  for  the  wave  number,  k.  The  growth 


rate  is  given  by  the  Im  k,  and 


Im  k 


1-1/2 


1 


(15) 


Note  that  the  growth  rate  is  always  higher  in  the  case  of  a  de- 

j|c 

celerated  beam  than  a  uniform  beam,  other  things  being  the  same. 

Finally,  it  may  be  added  that  by  making  measurements  at  three  dif¬ 
ferent  frequencies,  02^,  and  say,  the  three  unknown  quantities, 
o)p,  f3,  and  a,  can  be  determined  experimentally. 


This  can  be  seen  readily  when  one  calculates  an  "average"  value  of  k, 
(k>,  over  some  distance  z  =  L  (not  too  large),  say  where  (k)  ~  ~  k(o)  + 
k(L)j 
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Part  II  -  A  More  Rigorous  Approach  (H.  H.  C.  Chang) 


SUMMARY 


The  stability  of  a  non-relativistic  beam  of  charged  particles  in¬ 
jected  into  a  cold  plasma  is  studied  by  means  of  the  macroscopic  equations 
which  have  been  so  useful  in  analyzing  microwave  tubes  in  the  small  sig¬ 
nal  approximation.  In  this  approach,  the  simplifying  assumptions  are  made 
that  all  quantities  of  interest  vary  only  in  the  z  direction  and  that  they 
consist  of  a  d.c.  component  and  a  much  smaller  a.c.  component.  In  par¬ 
ticular,  1^,  the  d.c.  component  of  the  beam  current,  is  assumed  to  be 
uniform  and  infinite  in  extent  in  the  z  direction,  while  v  ,  the  d.c.  com- 
ponent  of  the  beam  velocity,  and  pQ,  the  d.c.  component  of  the  charge 
density,  vary  with  z  but  in  such  a  manner  that  Iq  =  P0VQ  is  constant.  To 
justify  the  neglect  of  the  boundary  effects  of  a  beam  of  finite  radius 
'a',  it  is  assumed  that  k^a  is  much  greater  than  unity,  where  k^  is  the 
transverse  propagation  vector.  For  low  beam  densities,  collisions  be¬ 
tween  the  beam  and  the  plasma  are  negligible.  The  plasma  is  assumed  to 
be  fully  ionized  and  only  collisions  between  the  electrons  and  the  ions 
are  taken  into  account  by  V,  a  velocity- independent  collision  frequency 
parameter.  Collisions  between  like  particles  and  thermal  effects  in  the 

cold  plasma  are  both  neglected.  In  Case  I,  it  is  assumed  that  v  = 

o 

V  exp(-az)  and  in  Case  II  that  v  =  V  [l  +  6  exp(-az)].  The  solution  of 
°  o  o 

the  dispersion  relation  is  discussed  in  some  detail  for  the  special  case 
of  wave  propagation  in  the  beam  direction.  The  general  case  when  the 
propagation  vector  makes  an  arbitrary  angle  with  the  beam  direction  is 
considerably  more  complicated,  and  with  certain  approximations  holding, 
the  dispersion  equation  is  shown  to  be  an  algebraic  equation  of  sixth 


10 


degree  in  cu,  the  wave  frequency.  Extension  to 
beam  is  discussed  briefly,  and  it  is  indicated 
upon  which  the  present  discussion  is  based,  is 
attainable  beam  currents. 


the  case  of  a  relativistic 
that  the  MHD  approximation, 
invalid  for  presently 
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In  this  paper,  the  stability  of  a  non- relativistic  beam  of  charged 
particles  injected  into  a  cold  plasma  is  studied.  The  beam  is  confined 
to  move  in  the  z  direction  only  and  is  infinite  in  extent  in  that  direc¬ 
tion.  It  is  uniform  in  the  xy  plane  and  its  radius  is  'a'.  Boundary 
effects  are  neglected  by  assuming  that  \j/a  is  much  less  than  unity, 
where  is  the  transverse  wavelength  [ 1 ] .  This  model  has  been  used  with 
considerable  success  by  tube  engineers  in  analyzing  microwave,  beam  type 
devices.  Invariably  a  sufficiently  strong  longitudinal  magnetic  field  is 
assumed  to  be  present  to  prevent  the  spreading  of  the  beam  due  to  space- 
charge  repulsion  of  the  electrons  of  the  beam.  To  keep  the  analysis 
tractable,  only  small  signal  phenomena  will  be  discussed.  This  means 
that  the  system  is  perturbed  very  slightly  from  an  equilibrium  state  and 
that  all  electrical  quantities  such  as  v(z),  the  beam  velocity,  I(z) ,  the 
beam  current  density,  p(z),  the  beam  charge  density,  and  the  electromag¬ 
netic  field  vectors  are  the  sum  of  a  d.c.  term  and  a  much  smaller  a.c. 
term.  Linearization  of  the  relevant  equations  is  justified  if  the  ratio 
(a.c.  term)/(d.c.  term)  is  always  negligible  compared  to  unity.  Col¬ 
lisions  between  the  beam  and  the  plasma  are  proportional  to  the  beam 
current,  and  for  the  low  beam  current  that  interests  us,  such  collisions 
may  be  neglected.  The  plasma  density  is  sufficiently  high  so  that 
electron-ion  collisions  within  the  plasma  must  be  taken  into  account  by 
V,  a  velocity-independent  collision  frequency  parameter.  Collisions 
between  life  particles  and  thermal  effects  in  the  plasma  are  neglected 
in  this  rough  calculation.  With  these  simplifying  assumptions,  we  are 
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studying  the  stability  against  small  displacements  of  a  beam  that  is 
penetrating  at  a  non-relativistic  velocity  v(z)  another  fully  ionized 
medium,  the  plasma.  This  model  patently  applies  to  microwave  beam  devices 
and  certain  astrophysical  and  thermonuclear  reactor  situations. 

I .  Non-Relativlstic  Beam 

We  begin  by  writing  the  pertinent  equations  governing  the  dynamics 
of  the  beam.  Let1 


I(z) 

=  Iq  +  I1(z)exp(jdbt)  , 

(i) 

p(z) 

=  pQ(z)  +  p1(z)exp(jo>t)  , 

(2) 

v(z) 

=  vq(z)  +  v1(z)exp(jtot)  , 

(3) 

with 

I  =  p  v  =  a  constant. 

0  0  0 

(4) 

In  the  above  equations,  the  subscript  zero  denotes  d.c.  quantities  and 
the  subscript  one  signifies  a.c.  quantities,  which  are  always  much 
smaller  than  the  corresponding  d.c.  quantities.  By  using  Maxwell's 
equations  and  the  Lorentz  force  law,  it  is  readily  shown  that  1^  satis¬ 
fies  the  following  equation:  [2] 


I."  +  (2jcu/v  +  3v'/v  ) I  *  +  (2ju>v'/v2  -  cu2/v2  -  el  /(emv3)  )I 
A  O  O  O  A  OO  OOOA 


=  -jewel  E/(mv  )  , 
o  o 


(5) 


where  e  is  the  permitivity  of  free  space  and  1^  =  dl^/dz. 


See  the  Glossary  for  the  definition  of  the  more  important 
in  this  report.  We  will  use  the  MKS  system  of  units. 

Eq.  (5)  seems  to  differ  from  Eq.  (8-35)  of  [3].  However, 
actually  identical.  If  in  (5)  we  replace  1^,  and  E  by 
and  E  -  jJ  / ea>  respectively,  we  will  obtain  flutter's  Eq. 

CiO  00 


symbols  used 


they  are 


J  ,  u 
cw  o 

(8-35) . 
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To  find  the  general  solution  to  (5)  let 


I^z)  =  Y(z)  exp(- jcoT)  =  Y(z)G(z)  =  YG 


where 


T  = 


dz/v 


Putting  (6)  and  (7)  into  (5)  we  get 


3  3 

Y"  +  3v'Y'/v  -  el  Y/(emv  )  =  -jtoel  E/(mv  G)  =  f(z) 

o  o  o  o  o  o 


(6) 

(7) 

(8) 


Let  us  first  find  the  homogeneous  solutions  of  (8).  Thus,  let 
y  (z)  and  y  (z)  be  two  linearly  independent  solutions  so  that  the 

i.  £• 

Wronskian,  W(z),  defined  by  (9)  is  non-zero. 


W(z)  =  y^y2 '  -  y2y^ '  . 


(9) 


The  particular  solution  of  (8)  is  then  [4] 

z  z 

yp(z)  =  -y f  (x)y2(x)/W(x)  dx  +  f  (x)y1(x)/W(x)  dx  (10) 


Case  I.  v  =  V  exp(-az) 
o  o 

Let  us  assume  that  due  to  collisions  and  other  dissipative  mechanisms 
the  beam  is  being  slowed  down  so  that 

v  =  V  exp(-az)  ,  0,  a  —  0.  (11) 

o  o 

From  (4)  it  follows  that 

Pq  =  exp(az)  ,  G(z)  =  exp(-ja>z/Vo)  .  (12) 
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Physically  this  assumption  is  meaningful  only  for  small  z  when  the 
approximation 

exp(±az)  =  1  ±  az  (13) 

is  a  good  one.  Upon  inserting  (11)  into  (8)  we  find  that  we  must  solve 
the  following  homogeneous  equation  for  y^(z)  and  y^z). 

Y"  -  3aY'  +  b  exp(3az)  Y  =  o,  (14) 

where  b  =  -el  /(emV3)  =  ne2/(emV2)  =  cr>2/V2  (15) 

o  o  o  p  o 

We  note  in  passing  that  if  a  =  0,  the  case  of  no  steady  electric  field 
or  energy  loss,  (14)  becomes 

Y "  +  u?  Y/V2  =  0  (16) 

P  o 

This  case  has  been  discussed  by  Smullin  in  Reference  2.  To  solve  (14) 
for  a  4  0,  set 

Y(z)  =  F(exp[£z] )  =  F(s),  s  =  exp(£z)  .  (17) 


Further,  if  £  =  3CZ/2,  (14)  becomes 


ds2 


1  dF  4bF 

s  ds  +  .2 

9a 


o,  a  6  o 


(18) 


In  Jahnke-Emde ,  p.  146,  the  solution  of  (18)  is  given  as 


F(s)  =  sZ  (a  s) ,  a  =  2to  /(3aV  )  ,  (19a) 

1  o  o  p  o 

where 

Z  (x)  =  AJ  (x)  +  BN  (x)  (19b) 

n  n  n 

is  the  cylindrical  function  and  J  (x)  and  N  (x)  are  the  Bessel  and 

n  n 
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Neumann  functions  of  order  n  respectively.  For  convenience  we  write  the 
solution  of  (14)  as 

Y(z)  =  uZ^(u)  ,  (20) 

where 

u  =  a  exp(3az/2)  .  (21) 

o 

It  follows  from  (19)  and  (20)  that  the  two  linearly  independent 
solutions  of  (14)  are 

yx(z)  =  uJ1(u)  and  y2(z)  =  ulVu*  (22) 

with  u  given  by  (21).  The  non-vanishing  Wronskian,  (9),  is 

W(z)  =  3au2/n  .  (23) 

With  vq  given  by  (11)  and  u  by  (21),  f(z)  of  (8)  can  be  written 

f(z)  =  -a  Eu2/G,  a  =  9  jcn2I  e/(4mu)2V  )  (24) 

i  1  o  p  o 

Putting  (22)  and  (24)  into  (10)  and  assuming  that  E  is  a  constant 
independent  of  z,  we  get 

na  Eu 

yp(z)  =  -55—  [J1(u)S2(z)  -  N1(u)S1(z)]  (25) 

where 

r  Z  (.1?)  d£  r  z  uN  (u)  dz 

S  (z)  =  I  - __  ,  S_(z)  =  — ~r—  (26) 

1  J  G(z)  2  G(z) 

It  is  easily  shown  that  v^(z)  is  obtained  from  I^(z)  in  the  fol¬ 
lowing  way:  [2] 
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(27) 


vx(z) 


V 

o 

I 

o 


(I^z)  -  j 


•  2 

JVoG  dY 

0)1  dz 
o 


V  I '  (z)/cd) 
o  1 


Upon  inserting  (11)  and  (20)  into  (27),  we  obtain  v  (z)  ,  the  homo- 
geneous  solution  of  the  differential  equation  which  v^(z)  satisfies. 

Vl  (z)  =  -  a  exp(az)G(z)[AJ  (u)  +  BN  (u)]  ,  (28) 

i.  z  o  o 

h 

where 

a  =  2j032/(3toI  a)  .  (29) 

&  p  o 

From  (6)  and  (20),  I  (z) ,  the  homogeneous  solution  of  (5),  is 

n 

(z)  =  uG(z)[AJ1(u)  +  BN1(u)]  .  (30) 


To  determine  the  constants  A  and  B,  we  define  z  and  z  by 

X  A 


z 

=  zi 

T  =  T 

u  =  ux 

,  N1(u1) 

=  o, 

(31) 

z 

=  Z2 

T  =  T2 

u  =  u2 

.  J  (u  ) 
o  z 

=  0  . 

(32) 

With  I  (z  )  and 
Ah 

v  (Zf) 

h 

known  from 

measurements , 

it  follows 

that 

A  =  l1  (z1)/[u1JJ(uJ)G(z1)]  ,  (33) 


B  =  "V1  (z2)exP(_az2)/^a2No(U2)G(Z2)-'  ‘  (34) 
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The  general  solution  of  (8)  Is 

Y(z)  =  u[AJ. (u)  +  BN.(u)]  +  y  (z)  ,  (35) 

lip 

where  v  (z)  is  given  by  (25)  and  u  =  a  exp(3az/2)  .  From  (6)  the 
'  P  ° 

general  solution  of  (5)  is 

Ix(z)  =  Y(z)G(z)  =  u[AJx(u)  +  BN1(u)}G(z)  +  g(z)E,  (36) 

where 

-  ita  uG(z) 

g(z)  =  - ^ -  [^(uJS^z)  -Jx(u)S2(z)]  ,  (37) 

and  S1(z)  and  S2(z)  are  defined  by  (26).  This  concludes  our  description 

of  the  beam,  and  we  will  now  describe  the  behavior  of  the  plasma. 

Let  the  cold  plasma  consist  of  electrons  of  mass  m  and  ions  of 

mass  M,  which  are  at  rest  in  equilibrium  with  densities  n  and  N  , 

o  o 

respectively.  In  the  absence  of  an  external  magnetic  field,  an  elec¬ 
tric  field,  E,  causes  a  current,  Ip,  to  flow  in  the  plasma,  where 

I  cm,2  2\  2  _ 

ja)IP  =  €l^~jv  +  “i!  1  =  ecDt  E  •  (38> 

cje>2  =  n  e2(m  +  M)/(emM)  o>2  =  (n  -  N  )2e2/[e(n  m  +  N  M)  ]  (39) 

e  o  i  o  o  o  o 

and  V  is  the  velocity-independent  collision  frequency  between  ions  and 

electrons . 

As  is  well  known,  E  satisfies  [6] 

_  _  , _  1~  1—  “ 

V  x[ V  x  E]  +  —  E  = - -  [1  +  I.  ]  .  (40) 

*  £t  P  D 

n  t-r\  r 
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In  (40),  I.  is  the  beam  current,  and  from  (36)  and  the  one-dimensional 
b 

model  we  are  using 

I  =  g(z)Ee  .  (41) 

0  z  z 

For  small  displacements,  the  equations  governing  the  dynamics  of 
the  beam  are  linear,  and  one  can  Fourier  analyze  the  motion  and  thus  re¬ 
gard  an  arbitrary  displacement  as  a  superposition  of  waves.  From  the 
value  of  the  displacement  and  of  its  rate  of  change  at  every  point  in 
space  at  a  given  time,  combined  with  the  dispersion  relation,  to  =  Oi(k) , 
one  can  then  calculate  the  subsequent  behavior  of  the  system. 

Let  us  seek  solutions  of  the  form  exp(j[uyt  -Tc  ■  r])  .  Formally 
we  may  replace  V  by  —  jTc  and  by  ja>.  With  these  substitutions  (40)  be¬ 
comes 

2  a?  “t2\  -  -  -  -  -  2 

k  -  --5  +  — )  E  -  (k  •  E)  k  +  jcoIb  /  (ec  )  =  0.  (42) 

c  c  / 

Orient  our  coordinate  axes  so  that  k  =  kxex  +  k  ez  *  Putting  (41) 

into  (42)  we  obtain  the  determinant  (43)  for  the  desired  dispersion  rela¬ 
tions. 

2  2  2  2  2 

ck-co+o)"  0  -ckk 

Z  t  X  z 

„  2(2  21  2  2  „  i 

0  clk+k  -oi+CD.  0 

1  x  z  I  t 

i 

2  2  2  2  2 
-  c  kxkz  0  c  kx  -  co  +  o>t  +  jcog(z)/e; 

(43) 

Eq.  (43)  can  be  solved  explicitly  for  the  case  of  wave  propagation 
in  the  direction  of  the  beam  (z)  axis,  when  kx  =  0.  In  this  special 
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case  we  find  that 


2  2  2,2 

to  =  o>t  t  c  k  ,  (44) 

is  a  double  root ,  and 

2  2 

co  =  oot  +  jcng(z)/e  (45) 

is  a  simple  root.  Prom  (15),  (24),  and  (37),  (45)  can  be  written 

co  =  cot[l  +  |  nauG(z)(N1(u)S1(z)  -  J^u^U)) ]~1/2  .  (46) 

We  are  interested  in  the  stability  of  the  beam-plasma  system  against 
small  disturbances  excited  at  t  =  0.  As  we  have  assumed  the  time  de¬ 
pendence  to  be  of  the  form  exp(jcot)  and  k  =  2n/\  is  assumed  to  be  real 
and  positive,  growing  waves  are  indicated  by  the  roots  to  =  to(k)  of  the 
dispersion  equation,  (43),  which  have  negative  imaginary  parts.  Eq.  (44) 
contains  no  reference  to  the  beam  parameters  and  relates  only  to  the 
properties  of  the  plasma.  In  particular,  it  contains  no  beam  convection 

term,  kv  =  kV  ,  and  if  V  =  0  it  has  real  roots  only.  If  n  =  N  and 
o  o  o  o 

2  2  2 

m/M  «  1,  co.  =0  and  to  =  co  =  n  e  /(em)  and  (44)  is  identical  with  the 
l  e  p  o 

dispersion  relationship  for  transverse  oscillations  in  a  plasma  7  . 

For  both  v/[co|  «  1  and  v/|co|  »  1,  it  can  be  verified  that  the  imaginary 
part  of  co  is  positive.  We  conclude  that  in  the  absence  of  the  beam,  the 
plasma  can  support  electromagnetic  waves.  This  is  physically  reasonable, 
for  the  dissipative  effects  of  collisions  in  the  plasma  can  only  lead  to 
the  damping  of  these  waves. 
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In  order  to  be  able  to  discuss  the  roots  of  (46)  we  must  obtain  an 
approximate  expression  for  the  denominator  of  the  right  hand  side  valid 
for  small  z.  Using  the  following  approximations, 


u  =  a  exp(3az/2)  =  a  (1  +  30fz/2) 
o  o 


J.(u)  =  J (a  )  +  [a  J  (a  )  -  J.(a  )] 

1  lo  2  ooo  lo 

„  2 
3aa  z 

uJ  (u)  =  a  J  (a  )  +  — —  J  (a  ) 

1  O  1  O  2  O  O 


(47) 


T  =  j  dz/vQ  =  i  exp(az)  dz  =  (eaz  -  1)  =  ~- 

Jq  0-0  O  O 


a  lengthy  but  elementary  calculation  shows  that  (46)  can  be  written 


approximately  as 
2 


tn2[l  + 


CD  .  CD 

— (  2(e"jPz  -  1)  +  jpz(3e“jpZ  -  2))]  =  co2  +  - — 


1805 


i  1  -  jV/05 


(48) 


2  2 

where  6  =  tn/V  and  to  =  ne  /em.  For  further  information  about  the  roots 

op 


of  (48),  we  require  the  following  approximation  to  be  valid: 


exp(-jpz)  =  1  -  jpz  (48') 

Substituting  (48')  into  (48),  the  dispersion  equation  becomes 

to3  -  (jv  +  jcu2  z/(18V  )  )co2  -  [to2  +  to2  +  zto2v/(18V  )]to+  juco2  =  0  (49) 

J  J  p  o  e  l  p  o  J  i 

In  the  absence  of  collisions  v  =  0,  and  (49)  has  solutions  to  =  0  and 

to  =  to  where  to  satisfies 
o  o 


2 

to 

o 


jzto2  to  /(18V  )  -  to2  -  to2 
p  o  o  e  i 


0  . 


(50) 
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2  2  2 

For  small  z,  0)  =  oi  +  m,  .  When  V/cjo  «  1,  we  can  write 

o  e  i  o 

co  =  coQ  (1  +  g/coo) ,  |/a>o  «  1  . 

Upon  substituting  this  into  (49)  and  using  (50)  we  learn  that 

2 

j\to 

co  =  co  +  - 5-  .  (51) 

o  „  & 

2co 

o 

The  last  result  means  that  for  wave  propagation  in  the  direction  of  the 

I  2  2 

beam  axis  and  when  v/coq  =  v/y(a>e  +  ux)  is  negligible  compared  to  unity, 
the  beam-plasma  system  is  stable  against  small  perturbations. 

In  a  fully  ionized  plasma,  three  types  of  interactions  between  the 
charged  particles  are  occurring:  electron-electron,  electron-ion,  and 
ion-ion  interactions.  We  have  assumed  for  simplicity  that  the  electron- 
ion  interactions  are  dominant  and  characterized  them  by  a  velocity  inde¬ 
pendent  collision-parameter  v.  The  collision  frequency  of  the  components 
of  a  plasma  is  a  complicated  subject  which  requires  careful  discussion. 

For  our  purposes,  we  follow  Delcroix's  book,  INTRODUCTION  TO  THE  THEORY 
OF  IONIZED  GASES,  where  on  page  111,  we  learn  that 

2 

V  =  4jtnw1(PQ  )avlog[h/pQ  ]  (52) 

av 

This  formula  shows  that  V  can  take  on  a  large  range  of  values  by  varying 
n,  the  number  density,  and/or  T,  the  temperature,  but  in  such  a  way  that 
we  are  always  dealing  with  a  non-relativistic ,  classical  plasma.  After 
specifying  cCp,  Vq  and  the  domain  of  z,  it  can  be  shown  that  by  varying 
V,  (49)  can  have  complex  roots  with  negative  imaginary  parts.  In  other 
words,  for  certain  values  of  v,  the  beam-plasma  system  is  unstable. 
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2  2  2 

In  the  general  case  when  k  =  k  +  k  ,  (43)  becomes 

x  x 

,  2  2  2,2xf/1  r>w  2  2W  2  2  2,2.  2, 2.  2  2V , 

(to  -  to  -  c  k  )  L  ( 1  +  P)(co  -  to  )(co  -  to  -  c  k  )  -  c  k  (to  -  to,)J  =  0, 

t  I  t  X  t 


or 


✓  2  2  2y2\t,,  TS\  4  ,  o  2  212  2  2,2lv  2 

(to  -o)  -ck)L(l  +  P)co  -  (2cjo  +  c  k  +  Puiy  +  c  k  ])o  + 

L  t  t  Z 


2 ,  2  22i  2  2  22 

cot(tot  +  c  k  )  ]  s  (to  -  cot  -  c  k  )F(to)  =  0 


(53) 


whe  re 


P(z)  =  P  =  -  rt3uG[N1(u)S1(z)  -  J1(u)S2(z)]  . 


(54) 


Note  that  when  k  =  k^,  (53)  yields  the  earlier  results,  (44)  and  (45), 
the  case  of  wave  propagation  in  the  direction  of  the  beam  axis,  as  is 
required.  Recalling  (12)  and  the  definition  of  P,  F(co)  =  0  is  a  tran¬ 
scendental  equation  in  to.  The  problem  is  to  determine  the  nature  of  its 
roots.  Using  (38),  the  explicit  form  of  F(co)  =  0  is 

F(to)  =  (1  +  P)tn6  -  2j V(  1  +  P)co5  -  ![1  +  P]v2  +  [2  +  p](co2  +  to2) 

»  e  l 

+  c2k2[l  +  P  cos2e]ioo4  +  j  V  ( [  2  +  P](to2  +  2co2)  +  2c2k2(l  +  Pcos20)ico3 

22  4  224222  2 

+  ico.V  (2  +  P)  +  co  +  2to  to.  +  to,  +  c  k  v  (1  +  P  cos  0) 
t  i  e  e  1  1 

2,2.  2  2.,  2  .  2  2  o  4  2,  2r  2  „  2nv 

+  c  k  (to  +  to.),  co  -  jV(2co  co.  +  2co.  +  c  k  [co  +  2ko.])to 
e  i  e  i  i  e  i 


2,  2  2  2  2  4 

-ckco.V  “  V  GO  =  0. 

1  1 


(55) 


2  2  2  2 

If  V  =  0,  co  =  co  =  co  +  co.  and  (55)  becomes 
t  t_  e  i 


F(co)  =  co*“[ (1  +  P)co**  -  ([2  +  P]to*  +  c2k2[l  +  P  cos‘s9])to“5  +  ay  (.uy  +  c^k2)]  =  0 

rO  to  to 

(56) 


2.2  2,  2V 
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Using  (6),  (19),  (21),  (26),  (47),  and  (54),  P(z)  can  be  approxi¬ 


mated  by 

2 

CD  . 

P(z)  i  [2(exP["j“z/v0i  "  1)  +  ir^(3  exP  [-Jcdz/V  ]-  2)]  (57) 

18cd  vo 

For  sufficiently  small  z  we  may  use  the  approximation 


exp(-jco/V  )  =  1  -  jcdz/V 

o  o 


(58) 


to  simplify  (57)  to  read 

.  2 

Jffi  z  „ 

p(2>  ;  -  Tisr  ■  -  I  •  ■>  -  *>p  (59> 

O 

Putting  (59)  into  (56)  we  learn  that  the  dispersion  equation  simplifies  to 


2r  4  3  2  2  2.  2  .  2  2,  2  2„,  2.2 

CD  I.CD  -  PCD  -  (2CDt  +  C  k  )(I3  +  (pCDt  +  pc  k  cos  0)CD  +  CDt  + 

o  o  o  o 


2,  2xl 

c  k  )J  = 


(60) 


According  to  a  celebrated  theorem  in  the  theory  of  equations  due  to  Euler 
and  Gauss,  "Every  polynomial  p(z)  of  positive  degree  with  complex  co¬ 
efficients  has  a  complex  root."  [8],  Eq.  (60)  fulfills  the  conditions  of 
this  theorem,  and  we  can  conclude  that  the  fourth  degree  polynomial  in  the 
square  bracket  has  a  complex  root.  The  solution  of  a  quartic  is  a 
straightforward,  albeit  laborious  process.  On  page  254  of  Reference  2, 

the  standard  recipe  for  solving  quartics  is  given.  It  is  patent  that  by 

2 

varying  k  =  2jt/\  and/or  p  =  jcd  z/18V  ,  (60)  can  be  made  to  possess  com- 

P  o 

plex  roots  co  =  cd^  ♦  jcd^  where  cd^  <  0.  In  varying  p,  however,  it  must  be 
remembered  that  z  must  be  sufficiently  small  so  that  the  series  expansions 
of  (47)  and  (58)  are  valid. 
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A  sensible  procedure  to  follow  in  our  search  for  the  nature  of  the 


roots  of  (56)  is  to  first  specify  a,  a>p,  and  Vq.  The  expansion 


f (a  +  z)  =  f (a  )  +  f'(a  )z 
o  o  o 


is  arbitrarily  defined  as  a  good  approximation  if 


(61a) 


f ’(a  )zi 
O  I 

f(a  ) 
o 


1000 


(61b) 


Let  z  ,  z  ,  z  ,  and  z  be  the  maximum  values  of  z  which  will  make  the  ex- 
12  3  4 

pansions  of  (47)  valid.  Let  zq  be  the  smallest  of  z^,  z z 3>  and  z^ 

In  (60)  we  then  allow  z  to  vary  between  z  =  0  and  z  =  zq.  Upon  speci¬ 
fying  0,  the  angle  between  k  and  the  z  axis  and  'a'  the  beam  radius,  we 
require  that 

k  sin  6  «  a  (62) 


and  investigate  if  there  are  any  complex  roots  of  (60),  satisfying  (63), 
whose  imaginary  parts  are  negative.  If  so,  the  system  is  unstable.  How¬ 
ever,  these  complex  roots  are  sensible  only  if  |cozo/Vo|  is  small  compared 
to  unity.  This  last  requirement  follows  from  our  use  of  the  approximate 
expression  exp(-joozo/Vo)  =  1  -  jcoz^/V^  to  obtain  (60). 

When  v  /  0  and  with  (59)  holding,  (55)  can  be  written 

to  -  (p  +  2jV)co  +  (2jvp  -  v  -2(coe  +  CD.)  -  c  k  )fl)  +  (plv  +  coe  +  CDi 


+  c"k2cos2e]  +  2jv[co2  +  2C02  +  2c2k2])co3  +  |(2co2  +  c2k2)v2 
e  l  i  i 

.  r  2  _  2  _  2,  2  2„1  .2  2.r  2,  2  2  2,,  2 

-jvplcoe  +  2coi  +  2c  k  cos  Gj  +  (co0  +  coi)lc  k  +  <jje  +  co^ J! CO 

2  2  „  4  2,  2r  2  „  2i  .2  .  2.  2  2Q. 

-jv(2co.co  +  2c o.  +  c  k  Leo  +  2co.J  -  jvco.p  -  jvc  k  p  cos  6)co 
i  e  i  e  i  l 

2  2  2  2  2.  //.q\ 

-  V  a>  (c  k  +  co^)  =  0  (63) 
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For  arbitrary  values  of  the  coefficients  of  cjo,  (63)  can  be  solved 
numerically  and  the  nature  of  the  roots  displayed.  The  procedure  out¬ 
lined  previously  for  the  case  v  =  0  can  be  followed  without  any  change. 
Although  this  procedure  is  simple  in  theory,  the  actual  numerical  work 
involved  can  be  quite  laborious  for  arbitrary  values  of  v/ct  and  will  not 
be  carried  out.  A  good  reference  on  this  subject  which  will  make  clear 
the  manifold  considerations  that  are  involved  is  Reference  3,  Chapters 

10  and  14.  This  concludes  our  discussion  of  the  case  where  v  -  V  exp(-az) 

o  o 

and  we  will  continue  by  discussing  a  different  dependence  of  v^  on  z. 

Case  II.  v  =  V  (1  +  ft  exp[-az]) 
o  o 

Assume  that 


vo  =  Vo(l  +  p  exp[-az]),  0  <  ft  «  1,  a  >  0  .  (64) 

As  we  require  that  I  =  p  v  =  constant  to  terms  of  order  ft,  we  must 
o  o  o  K 

put 


PQ  =  PQ(1  -  ft  exp[  -az] )  . 


(65) 


I^z),  the  a.c.  current,  satisfies  (5).  If  we  put 


r  /  \  Y(T)  ,  .  GY(T) 

1 1  ( z)  =  exp  (-ju>T)  =  y 


(66) 


where 


/  d  z 

T  =  /  —  ,  (5)  becomes 

o 


d2Y 

dT2 


1  <2v 

_1 _ o 

V 

o  dT 


el  -jecDl  E 

o'  o 

- 2  I  Y  =  - 

mv  !  mG 


fx(T) 


(67) 


Neglecting  ft  compared  to  unity  and  with  (13)  holding,  (67)  can  be  approx¬ 
imated  by 
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fx(T)  , 


(68) 


where 


-  (a_exp(-aV  T)  +  o>  2)y  = 
dT  op 


2  2 

a3  =  pa  VQexp  p 


_  2  2 
s  pa  v 

^  o 


CD 


-eI/(emV  ) 
o  o 


As  before,  let  us  first  find  the  homogeneous  solutions  of  (68).  To 
end,  let 

Y(T)  =  F(exp[-£T])  =  F(s),  s  =  exp[-£T]  ,  £  =  aVQ 


It  then  turns  out  that  F(s)  satisfies 


2  d2F  dF  1  .  2. 

S  ~7.  +  S  di  “  12  (a3S  +  V*  =  °- 

as  £ 


According  to  page  146  of  Jahnke-Emde,  solution  of  (71)  is 


F(s)  =  AJn(2VF7  )  +  BNn(2YpT),  n  =  2cop/£ 

The  Wronskian,  W(T)  as  defined  by  (9)  of  y^(T)  =  Jn(2-\/ps  )  and 
y2(T)  =  Nn(2\/ps  )  is 


W(T)  = 


.  L 


The  particular  solution  of  (68)  is,  according  to  (10), 

T 


yp(T)  = 


2  „ 

-jjtuxd  eE 

P 

a 


N  (9)  j  J  (0)exp( jcoT)  dT  - 
n  Jo  n-  -  - 


J  (9)  /  N  (6)  exp(jCDT) 
n  J  n  —  — 


-jnaiD  eE 


dT  s 


a 


W(9) 


where  9  =  2"\/ps 


(69) 

this 


(70) 


(71) 


(72) 


(73) 


(74) 
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The  complete  solution  of  (5) ,  with  (64)  holding,  is 


I^z) 


~  LAJn(9)  +  BNn(0)j  ♦  gl(z)E, 


(75) 

where 

-jnrnc  eG 

gi(z)  =  ^7^-w(e) 

o 

The  homogeneous  solution  to  the  equation  which  v^(z)  satisfies 

will  be  designated  by  v  (z),  as  before.  It  follows  from  (27)  and  (75) 

xh 


that 

2 

jav  GG 

vi  <z>  =  “5^7 -  [AJ'(8)  +  BN' (9)  ]  ,  (76) 

l,  zcX)i  n  n 

h  o 

where  J'(9)  =  dj  (9)/d9  .  I,  (z) ,  the  homogeneous  solution  of  (5),  is 

n  n  1 


Ix  (z)  =  7-  [AJn(0)  +  BN  (9)]  .  (77) 

n  o' 

To  determine  the  constants  A  and  B,  we  define  z  and  z  by 

X  £i 

z  =  zx,  T  =  Tx,  9  =  01,  VV  =  ° 

z  =  z2,  T  =  T2,  9  =  02,  J*(02)  =  0 


With  I  (z  ) 
h 


and  (z^)  known  from  measurements,  it  follows  that 
h 


A  =  I  (z  )V  /[G(T  )J  (9  )] 
1.  1  o  1  n  1 

h 

B  =  2cdI  v  (z)/[0  V2G(T  )] 

0  1,  i  O  i 
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The  interaction  between  the  beam  and  the  plasma  is  accounted  for 
exactly  as  in  the  previous  case.  Thus,  the  analysis  following  Eq.  (37) 
is  applicable  with  slight  changes  in  notation.  Eq.  (41)  becomes 


h  =  gl(2)Vz  •  <79> 

and  (43)  still  holds  with  g^z)  replacing  g(z).  Eq.  (44)  is  unchanged 
and  (46)  becomes 


[i  +  jtoj2  G(z)w(e)/(av  )]1/2 

P  O 


(80) 


In  the  general  case,  when  k 


2  2 

k^  +  kz  ,  we  get  an  equation  of  the 


form  given  by  (55)  except  here  we  must  replace  P(z)  by  P^(z)  where 


P.(z)  =  net  G(z)W(0)/(QV  ) 

1  p  o 


(81) 


For  sufficiently  small  z,  a  and  p  we  can  obtain  an  approximate  ex¬ 
pression  for  W(0)  as  defined  by  (74).  Thus  by  means  of  series  expansions 

similar  to  (47)  we  learn  that  W(0)  =  -QV  T2/(2n)  and 

o 


-  U) 


P^z)  = 


P  (1  +  2|3)  ( 1  +  [1  +  p]  4^  >*2 


(82) 


2V 


When  jwz(l  +  P)/Vq  is  much  less  than  unity,  we  approximate  P^z)  by 


-  (1  +  2P)z2w  2 

P,(z)  =  - 5 - 2_  (83) 

1  2  V  2 

o 

Putting  P - >  P^(z)  in  (55)  we  get  the  sixth  degree  dispersion  equation 

which  must  be  solved  in  order  to  predict  if  the  beam-plasma  system  is 
stable  or  not.  As  in  Case  I,  this  sixth  degree  equation  can  be  solved 
numerically  after  specifying  the  parameters  which  relate  to  the  properties 
of  the  beam  and  the  plasma. 
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II.  Relativistic  Beam 


The  rest  mass  of  an  electron  is  about  one-half  MEV.  The  kinetic 
energy  of  an  electron  in  rest  mass  units  is  7  -  1,  and  when  y  -  1  is 
much  greater  than  unity,  the  previous  considerations  must  be  modified 
because  the  neglected  relativistic  corrections  are  then  important.  For 
a  relativistic  beam,  examination  of  the  equations  upon  which  our  crude 
one-dimensional  model  is  based  shows  that  the  linearized  equation  of  con¬ 
tinuity  does  not  change  and  there  is  no  change  in  Poisson's  equation 

unless  a>/k  =  c.  [9].  To  the  first  approximation,  all  we  need  to  do  to 

3 

include  relativistic  effects  is -to  replace  m  by  1 17  as  far  as  the  beam 
is  concerned.  A  more  sophisticated  and  accurate  treatment  based  on 
Boltzmann's  equation  is  given  in  Reference  5. 

In  truth  for  extremely  energetic  beams,  say  when  y  -  1  is  much 
greater  than  100,  the  situation  is  much  more  complicated  than  just  indi¬ 
cated.  For  one  thing,  as  Budker  has  shown,  for  highly  relativistic 
beams,  radiation  is  an  important  dissipative  mechanism  [lo].  Moreover, 
if  the  beam  energy  exceeds  the  relevant  threshold  energies,  other  ele¬ 
mentary  particles  such  as  mesons  could  be  produced  copiously.  Additional¬ 
ly  this  entire  discussion  is  essentially  a  magnetohydrodynamic  one  and 
that, as  Finkelstein  and  Sturrock  [ 11 ]  have  shown  for  a  relativistic  beam, 
the  MHD  approximation  is  valid  provided  the  number  of  electrons  per 
classical  electron  radius,  measured  in  the  rest  frame  of  the  rapidly 
moving  electrons,  is  large  compared  with  unity.  In  practice,  this  last 
requirement  is  equivalent  to  the  statement  that  the  current  of  the  beam 
must  exceed  17,000  7  amps  and  this  current  is  much  greater  than  currents 
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so  far  contemplated  in  relativistic-stream  experiments. 

In  Reference  5,  the  stability  of  a  uniform,  relativistic  beam  of 

particles  injected  into  a  plasma  is  studied  using  the  Boltzmann  transport 

equation.  The  most  important  results  which  were  obtained  are:  (1)  The 

beam  is  relativistic  so  that  in  the  equations  of  motion  the  relativistic 
3 

masses  and  y  m^  appear.  This  makes  the  beam  "stiff"  so  that  magnetic 
self-forces,  which  for  non-relativistic  pinches  lead  to  kinking,  are 
negligible  compared  with  electrostatic  forces.  (2)  The  velocity  spread 
in  the  beam  is  allowed  to  be  appreciable  so  that  for  short  enough  wave¬ 
lengths  the  velocity  structure  of  the  beam  is  important.  The  velocity 
spread  in  the  plasma  penetrated  is  small.  (3)  As  in  this  report,  only 
electron- ion  collisions  are  considered  important  and  are  taken  care  of  via 
a  velocity  independent  collision  frequency  parameter.  In  this  important 
paper,  the  effects  of  boundary  conditions  in  a  finite  beam  and  the  non- 
uniform  structure  present  when  velocity  gradients  are  significant  are 
assumed  to  be  small  and  consequently  negligible. 

In  Reference  12,  two-stream  instability  in  finite,  relativistic 
streams  is  investigated.  The  electrostatic  mode  is  singled  out  for  dis¬ 
cussion,  and  a  detailed  study  is  made  of  effects  associated  with  the 
finite  size  and  non-uniformity  of  the  beam.  An  eigenvalue  equation  for 
Ez,  the  longitudinal  electrical  field,  rs  obtained  and  for  a  weak-beam 
it  is  shown  that  two-stream  instability  can  exist  only  if  t]  is  real  and 
T|  S  !•  Here  t|  is  a  function  of  co,  v,  and  the  plasma  frequency  of  the 
beam  and  is  defined  by  Eq.  (29)  of  this  reference. 
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Finally  we  note  that  in  some  ways  a  study  of  the  instabilities  of  a 
relativistic  beam  is  simpler  than  that  of  a  non-relat ivist ic  beam.  In 
using  a  one-dimensional  model  to  describe  a  non-relativistic  beam  it  is 
imperative  that  a  large  longitudinal  d.c.  magnetic  field  be  present  to 
inhibit  transverse  motion  of  the  electrons  and  thus  keep  the  beam  from 
spreading.  An  extreme  relativistic  beam  is  pinched  by  its  own  magnetic 
field  and  when  the  Bennett  condition  [11]  is  satisfied,  self-magnetic 
effects  are  sufficient  to  hold  the  beam  together  against  its  random 
transverse  kinetic  energy  and  multiple  scattering  effects  without  an  ex¬ 
ternally  applied  longitudinal  magnetic  field.  Obviously  the  complexity 
of  the  analysis  is  strongly  dependent  on  the  parametric  choices  made  as 
to  whether  the  skin  depth  is  large  or  small  compared  to  the  beam  radius 
and  the  magnitude  of  the  ratio  |v/u>|.  The  various  instabilities  are  all 
means  of  transferring  energy  from  the  beam  to  the  plasma  substratum,  and 
they  divide  naturally  into  those  which  do  and  which  do  not  depend  on 
edge  effects — whether  the  beam  has  sharp  edges  or  diffuse  edges.  A  good 
reference  on  the  stability  of  a  relativistic  beam  is  Reference  13. 
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GLOSSARY 


a 

a 

o 


b 


c 

e 


e 

z 

E 

F 

g(z)  =  g 
G(z)  =  G 
I 


I^z) 


=  I, 


Jn(z) 

k 

k,, 


beam  radius 

constant  defined  by  Eq.  (19a) 
constant  defined  by  Eq.  (24) 
constant  defined  by  Eq.  (29) 
constant  defined  by  Eq.  (69) 
b  =  Oip/V^  ,  Eq.  (15) 
velocity  of  light 

magnitude  of  the  charge  of  the  electron, 

-19 

e  =  +  1.6  x  10  coulomb 

unit  vector  in  the  z  direction,  the  beam  direction 
intensity  of  the  electric  field  vector 
function  of  z  defined  by  Eq.  (17) 
function  of  z  defined  by  Eq.  (37) 

G  =  exp(-jcuT) 

beam  current,  I(z)  =  Iq  +  I^z)  exp(jcot) 
d.c.  component  of  beam  current 
a.c.  component  of  beam  current 
1^  -  dl  /dz 

2 

imaginary  unit,  j  =  -1 

Bessel  function  of  order  n 

propagation  vector 

z  component  of  k,  kn  =  k  cos  9 

normal  component  of  k,  kx  =  k  sin  0 

electron  mass 

ion  mass 
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n 

o 

number  density  of  electrons  of  the  plasma 

N 

o 

number  density  of  ions  of  the  plasma 

N„(z) 

Neumann  function  of  order  n 

P 

quantity  defined  by  Eq.  (59) 

P(z)  =  p 

quantity  defined  by  Eq.  (54) 

s 

variable  defined  by  Eq.  (17),  s  =  exp(£z) .  Also 
s  =  exp(-jT) 

S  ,  S 

1’  2 

quantities  defined  by  Eq.  (26) 

t 

time 

r  Z 

T 

T  =  /  dz/v 

J 

u 

o 

variable  defined  by  Eq.  (21),  u  =  aoexp(3az/2) 

_  _ 

v(z)  =  V 

velocity  of  the  beam,  v  =  v  e^  , 
v  =  vq  +  v  (z)  exp(jo)t) 

V  (z)  = 

O  0 

d.c.  component  of  v 

vl(z)  =  v1 

a.c.  component  of  v 

V 

o 

magnitude  of  vq  in  Case  I  at  z  =  0;  quantity  defined 

in  Eq.  (64)  related  to  the  magnitude  of  v 

o 

W(z) 

Wronskian,  Eq.  (9) 

z 

the  independent  variable 

Z 

n 

the  cylindrical  function,  Eq.  (19b) 

a 

exponential  decrement  factor,  Eq.  (11)  and  Eq.  (64) 

P 

in  Case  I,  g  =  co/V  .  (3  is  also  used  as  in  Eq.  (64) 

;  va 1-1/2  ° 

7 

Z  -  I1  -  -2 

'  c  1 

e 

permitivity  of  free  space 

5 

II 

Q 

< 

o 

e 

angle  between  k  and  e 

z 
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wavelength,  k 


2n/\ 


i 


\ 


p(z) 


03 


03. 

1 

03 

P 

03 

t 

03 

o 


=  \  sin  © 

\u  =  \  cos  0 

collision  frequency  between  electrons  and  ions,  Eq.  (52) 
p  electric  charge  density,  p  =  p^  +  p^ 
d.c.  component  of  p 
wave  frequency 
quantity  defined  by  Eq.  (39) 
quantity  defined  by  Eq.  (39) 
plasma  frequency,  Eq.  (15) 

quantity  defined  by  Eq.  (39) 

2.2  2 

03  =  03  +  03. 

o  e  i 
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